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d ' Traditionalcomputers work with finite numbers. Situations where the us- 

age of infinite or infinitesimal quantities is required are studied mainly the- 
oretically. In this paper, a recently introduced computational methodology 
(that is not related to the non-standard analysis) is used to work with finite, 
infinite, and infinitesimal numbers numerically. This can be done on a new 
»^ , kind of a computer - the Infinity Computer - able to work with all these types 

\^ ' of numbers. The new computational tools both give possibilities to execute 

computations of a new type and open new horizons for creating new mathe- 
matical models where a computational usage of infinite and/or infinitesimal 
Cn . numbers can be useful. A number of numerical examples showing the poten- 

^D ■ tial of the new approach and dealing with divergent series, limits, probability 

^^ , theory, linear algebra, and calculation of volumes of objects consisting of 

parts of different dimensions are given. 
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1 Introduction 

The point of view on infinity accepted nowadays takes its origins from the famous 
ideas of Georg Cantor (see ||2l). Different generalizations of the traditional arith- 
metic for finite numbers to the case of infinite and infinitesimal numbers have been 
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proposed in literature (see, e.g., lH] |2l |3] HI |3 and references given therein). How- 
ever, these generalizations are quite different with respect to the finite arithmetic 
we are used to deal with. Moreover, very often they leave undetemiined many op- 
erations where infinite numbers take part (e.g., °° — °°,^, etc.) or use representation 
of infinite numbers based on infinite sequences of finite numbers. 

In spite of these crucial difficulties and due to the enormous importance of the 
concept of infinite and infinitesimal in Science, people try to introduce these no- 
tions in their work with computers. Thus, the IEEE Standard for Binary Floating- 
Point Arithmetic (IEEE 754) being the most widely-used standard for floating- 
point computation defines formats for representing special values for positive and 
negative infinities and NaN (Not a Number) (see also incorporation of these notions 
in the intei-val analysis implementations e.g., in [ilSil ). The IEEE infinity values can 
be the result of arithmetic overflow, division by zero, or other exceptional opera- 
tions. In turn, NaN is a value or symbol that can be produced as the result of a 
number of operations including that involving zero, NaN itself, and infinities. 

Recently, a new applied point of view on infinite and infinitesimal numbers 
has been introduced in ISKIlJiniOini. With respect to the IEEE 754 standard, 
the new approach significantly extends the variety of operations that can be done 
with infinity. It gives a possibility to work with various infinite and infinitesimal 
quantities numerically by using a new kind of a computer - the Infinity Computer 
- introduced in ||9j \T0\ [Til. A number of applications related to the usage of the 
new approach for studying fractals (being one of the main scientific interests of the 
author (see, e.g., |[T6l[T7l ) has been discovered (see |[T3]fT4l '). 

The new approach is not related to non-standard analysis ideas from [Tl and 
does not use Cantor's ideas (see |13) either. The Infinity Computer works with 
infinite and infinitesimal numbers numerically using the following methodological 
principles having a strong applied character (see survey [fTSl for a detailed discus- 
sion on the new approach). 

Postulate 1. Existence of infinite and infinitesimal objects is postulated but it 
is also accepted that human beings and machines are able to execute only a finite 
number of operations. 

Postulate 2. It is not discussed what are the mathematical objects we deal 
with; we just construct more powerful tools that allow us to improve our capacities 
to observe and to describe properties of mathematical objects. 

Postulate 3. The principle formulated by Ancient Greeks "Hie part is less than 
the whole ' is applied to all numbers (finite, infinite, and infinitesimal) and to all 
sets and processes (finite and infinite). 

Due to this declared applied statement, such traditional concepts as bijection, 
numerable and continuum sets, cardinal and ordinal numbers are not applied when 
one works with the Infinity Computer because they belong to Cantor's approach 
having significantly more theoretical character and based on different assumptions. 
However, the methodology used by the Infinity Computer does not contradict Can- 
tor. In contrast, it evolves his deep ideas regarding existence of different infinite 
numbers in a more practical way. 



The accepted applied methodology means in particular (see Postulate 1) that we 
shall never be able to give a complete description of infinite processes and sets due 
to our finite capabilities. Acceptance of Postulate 1 means also that we understand 
that we aie able to write down only a finite number of symbols to express numbers 
in any numeral systenu. 

Postulate 2 states that the philosophical triad - researcher, object of investi- 
gation, and tools used to observe the object - existing in such natural sciences as 
Physics and Chemistry, exists in Mathematics, too. In natural sciences, the instru- 
ment used to observe the object limits and influences the results of observations. 
The same happens in Mathematics where numeral systems used to express num- 
bers are among the instruments of observations used by mathematicians. Usage 
of a powerful numeral system gives a possibility to obtain more precise results in 
Mathematics in the same way as usage of a good microscope gives a possibility 
to obtain more precise observations in Physics. However, due to Postulate 1, the 
capabilities of the tools will be always limited. 

Particularly, this means that from this applied point of view, axiomatic systems 
do not define mathematical objects but just deteiTnine formal rules for operating 
with certain numerals reflecting some properties of the studied mathematical ob- 
jects. For example, axioms for real numbers are considered together with a particu- 
lar numeral system S used to write down numerals and are viewed as practical rules 
(associative and commutative properties of multiplication and addition, distributive 
property of multiplication over addition, etc.) describing possible operations with 
the numerals. The completeness property is interpreted as a possibility to extend S 
with additional symbols (e.g., e, %, \/2, etc.) taking care of the fact that the results 
of computations with these symbols agree with facts observed in practice. As a 
rule, assertions regarding numbers that cannot be expressed in a numeral system 
are avoided (e.g., it is not supposed that real numbers form a field). 

In this paper, the methodology from ||8] [TT] [13] [14] [15] is used to describe how 
the Infinity Computer can be applied for solving new and old (but with higher pre- 
cision) computational problems. Representation of infinite and infinitesimal num- 
bers at the Infinity Computer and operations with them are described in Section |2] 
Then, Sections [3] and |4] present results dealing with applications related to hnear 
algebra and calculating divergent series, limits, volumes, and probabilities. 

2 Representation of infinite and infinitesimal numbers at 
tlie Infinity Computer 

In |[8][ll][T3][T4][T5l, a new powerful numeral system has been developed to ex- 
press finite, infinite, and infinitesimal numbers in a unique framework by a finite 
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number of symbols. The main idea of the new approach consists of measuring infi- 
nite and infinitesimal quantities by different (infinite, finite, and infinitesimal) units 
of measure. This section gives just a brief tour to the representation of infinite and 
infinitesimal numbers at the Infinity Computer and describes how operations with 
them can be executed. It allows us to introduce the necessary notions and designa- 
tions. In order to have a comprehensive introduction to the new methodology, we 
invite the reader to have a look at the recent survey [|T5l downloadable from ifTOl or 
at the book lH (written in a popular manner) before approaching Sections [3] and IH 

2.1 A new infinite numeral and a positional numeral system with 
the infinite radix 

A new infinite unit of measure has been introduced for this purpose as the number 
of elements of the set N of natural numbers. It is expressed by a new numeral (D 
called grossone. It is necessary to emphasize immediately that the infinite number 
® is not either Cantor's No or CO. In particular, ® has both cardinal and ordinal 
properties as usual finite natural numbers. Formally, grossone is introduced as 
a new number by describing its properties postulated by the Infinite Unit Axiom 
(lUA) (see 18] [TOl [121 [IS [13). This axiom is added to axioms for real numbers 
(viewed in sense of Postulates 1-3) similarly to addition of the axiom determining 
zero to axioms of natural numbers when integer numbers are introduced. 

One of the important differences of the new approach with respect to the non- 
standai^d analysis consists of the fact that (D G N because grossone has been intro- 
duced as the quantity of natural numbers (similarly, the number 3 being the number 
of elements of the set {1,2,3} is the largest element in this set). The new numeral 
® allows one to write down the set, N, of natural numbers in the form 

N = {1,2,3, ... ®-3, ®-2, ®-l, ®] (1) 

where the numerals 

... (l)-3, ®-2, ®-l, ® (2) 

indicate infinite natural numbers. 

It is important to emphasize that in the new approach the set ([T|) is the same set 
of natural numbers 

N = {1,2,3, ... } (3) 

we are used to deal with and infinite numbers (O also belong to N. Both records, 
([B and ([3]), are correct and do not contradict each other. They just use two different 
numeral systems to express N. Traditional numeral systems do not allow us to 
see infinite natural numbers that we can observe now thanks to ®. Similarly, a 
primitive tribe of Piraha (see ||4]) living in Amazonia that uses a very weak numeral 
system for counting (one, two, many) is not able to see finite natural numbers 
greater than 2. In spite of this fact, these numbers (e.g., 3 and 4) belong to N and 
are visible if one uses a more powerful numeral system. Thus, we have the same 



object of observation - the set N - that can be observed by different instruments - 
numeral systems - with different accuracies (see Postulate 2). 

It is worthy to notice that the introduction of ® defines the set of extended 
natural numbers indicated as N and including N as a proper subset 

N = {l,2,...,(l)-l,®,®+l,...,®2-l,®^®2 + l,...}. (4) 

Due to Postulates 1 and 2, the new numeral system cannot give answers to all 
questions regarding infinite sets. What can we say, for instance, about the number 
of elements of the set N? The introduced numeral system based on ® is too weak 
to give an answer to this question. It is necessary to introduce in a way a more 
powerful numeral system by defining new numerals (for instance, ®, ®, etc). 

Inasmuch as it has been postulated that grossone is a number, associative and 
commutative properties of multiplication and addition, distributive property of 
multiplication over addition, existence of inverse elements with respect to addi- 
tion and multiplication hold for grossone as for finite numbers. Particularly, this 
means that the following relations hold for grossone, as for any other number 

0-® = ®.0 = 0, ®-® = 0, ^ = 1, ®'' = 1, 1® = 1, 0® = 0. (5) 

To express infinite and infinitesimal numbers at the Infinity Computer, records 
similai" to traditional positional numeral systems can be used (see ll8l [T0l[T2l[T3l ). 
In order to construct a number C in the new numeral positional system with the 
radix ®, we subdivide C into groups corresponding to powers of ®: 

C = Cp,„®"'" + . . . + Cp, ®P' + Cp„®"» + Cp_, ®P-' + ...+ Cp_^®P-' . (6) 

Then, the record 

C = cp,„®^"' . . . cp, ®'^'cp„®^<'cp_, ®^-' . . . cp_^®>'-' (7) 

represents the number C, where c, are called grossdigits and are expressed by tra- 
ditional numeral systems used to represent finite numbers (e.g., floating point nu- 
merals). Grossdigits can be both positive and negative. They show how many 
corresponding units should be added or subtracted in order to foim the number C. 
Grossdigits can be expressed by several symbols. 

Numbers p, in (jT]) called grosspowers can be finite, infinite, and infinitesimal. 
They are sorted in decreasing order with po = 0- 

Pm > Pm-l > ■■■> Pl> Pq> P-l> ■■■P-(k-\) > P-k- 

In the record (|7]), we write ®^'' explicitly because in the new numeral positional 
system the number / in general is not equal to the grosspower pi. 

Finite numbers in this new numeral system are represented by numerals having 
only the grosspower po = 0. In fact, if we have a number C such that m = k = 
in representation (|7]i, then due to ^, we have C = co®° = cq. Thus, the number 



C in this case does not contain grossone and is equal to the grossdigit cq being a 
conventional finite number expressed in a traditional finite numeral system. 

Infinitesimal numbers are represented by numerals C having only negative fi- 
nite or infinite grosspowers, e.g., 3.48(1)^^^'^'26.4(1)^'^^®. The simplest infinitesi- 
mal number is tjn = ® being the inverse element with respect to multiplication 
for®: 

®-K® = ®-®-^ = 1. (8) 

Note that all infinitesimals are not equal to zero. Particularly, 4^ > because it is 
a result of division of two positive numbers. 

Infinite numbers are expressed by numerals having at least one positive finite 
or infinite grosspower. Thus, they have infinite parts and can also have a finite part 
and infinitesimal ones. For instance, the number 

12.4®^'^-2'®(-20.64)®'^0.8®°0.71®"^32.1®"^-^® 

has two infinite parts, 12.4®-''^-^^® and -20.64®^^ one finite part, 0.8®° = 0.8, 
and two infinitesimal parts, 0.71®"^ and 32.1®"^'^®. 

2.2 Arithmetical operations executed by the Infinity Computer 

A working softwai^e simulator of the Infinity Computer has been implemented (see 
||9j[T0l[TTl). It works with infinite, finite, and infinitesimal numbers numerically, 
(not symbolically) and executes the arithmetical operations as follows. 

Let us consider the operation of addition {subtraction is a direct consequence 
of addition and is thus omitted) of two given infinite numbers A and B, where 

KM L 

A = £fl,,®'^', B=£ft,„^®"^ C = £q,®'', (9) 

i=\ 7=1 i=\ 

and the result C = A + B is constructed by including in it all items a,t,® ' from A 
such that ki ^mj,l < j < M, and all items bm ®"^' from B such that mj ^ki,l < 
i < K. If in A and B there are items such that ki = m,, for some / and j, then this 
grosspower ki is included in C with the grossdigit b^. +aki, i.e., as {b^. +ak.)® '. 

The operation of multiplication of two numbers A and B in the form ^ returns, 
as the result, the infinite number C constructed as follows: 

M K 

C=LC;, C,• = ^„^®'"^•A = £ai,^,^.®'^'+'"^ \<j<M. (10) 

7=1 i=\ 

In the operation of division of a number C by a number B from Q, we obtain 
a result A and a reminder R (that can be also equal to zero), i.e., C = AB + R. The 
number A is constructed as follows. The first grossdigit ajt^ and the corresponding 
maximal exponent kx are established from the equalities 

akK = ci^lhrnu 1 ^K = h- mu- (1 1) 



Then the first partial reminder Ri is calculated as 

Ri=C-ak^®'"' -B. (12) 

If Ri 7^ then the number C is substituted by R\ and the process is repeated by a 
complete analogy. The process stops when a partial reminder equal to zero is found 
(this means that the final reminder R = 0) or when a required accuracy of the result 
is reached. 

Example 2.1. We consider two infinite numbers A and B, where 

A = 304.21®^''-^*(-7.1)®'241.2®*', B = 6.23®^13.1®*'l5®"^-2*. 

Their sum C is calculated as follows: 

C=A+S = 304.21®^^-^* + (-7.1)®^^ + 41.2®° + 6.23®V 13.1®° + 15®"^-^'" = 
304.21®^^-^* - 7.1®^2 _^ 6.23®^ + 54.3®" + 15®"^-^® = 

304.21®'^-*^®(-7.1)®^26.23®^54.3®°15®"^-2®. □ 

More examples illustrating the work of the Infinity Computer can be found in 

iiiiiiiniiniiia. 

3 Examples of situations where the Infinity Computer 
executes operations that traditionally required a human 
intervention 

In this section, we describe a number of computational tools provided by the new 
methodology and the Infinity Computer. It becomes possible in several occasions 
to automatize the process of the solving of computational problems avoiding an 
inten^uption of the work of computer procedures and the necessity of a human in- 
tervention required when one works with traditional computers. For instance, when 
one meets a necessity to work with divergent series or, even worse, their difference, 
traditional computers are not able to execute these operations automatically and a 
human intervention is required in order to avoid these difficulties. In this situation 
and in other examples below, it is shown how the work usually done by humans 
can be formalized following Postulates 1-3 and passed to the Infinity Computer. 

It is necessary to emphasize that the examples described in this section are 
related to numerical computations at the Infinity Computer. No symbolic compu- 
tations are required to work with infinite and infinitesimal numbers when one uses 
the Infinity Computer. 



3.1 Calculating sums with an infinite number of items 

The new approach allows one to use the Infinity Computer for the calculation of 
sums with an infinite number of items. The term 'series' is not used here because, 
due to Postulate 3, it is required to indicate explicitly the number of items (finite 
or infinite) in any sum. Naturally, it is necessary that the number of items and 
the result of the considered sum ai^e expressible in the numeral system used for 
calculations. 

Let us illustrate the new possibilities by considering two traditional infinite 
series Si = 1 + 1 + 1 + ... and ^2 = 30 + 30 + 30 + . . . The traditional analysis gives 
us a very poor answer that both of them diverge to infinity and, therefore, the results 
cannot be calculated and represented at the traditional computers. Such operations 
as ^2 — Si or 1^ are not defined and humans should return to the original physical 
problem in order to understand whether there exist answers to these questions. 

Now, when we are able to express not only different finite numbers but also 
different infinite numbers, it is necessaiy to indicate explicitly the number of items 
in the sums Si and S2 and it is not important whether it is finite or infinite. Due to 
Postulate 3, by changing the number of items in the sums, we change the respective 
results, too. 

Suppose that the sum Si has k items and the sum S2 has n items. Then 

Si (A:) = 1 + 1 + 1 + ... + !, S2(?i)=30 + 30 + 30 + ... + 30. 



and it follows Si {k) = k and S2(«) = 30«. If the results Si {k) = k and S2(«) = 30« 
are expressible in the chosen numeral system, then indeterminate forms disappear 
and the expressions S2(^) — Si {n) and ^44 can be easily calculated using the In- 
finity Computer. 

If, for instance, k = n = 5® then we obtain Si (5®) = 5®, S2(5®) = 150® and 
S2(5®)-Si(5®) = 145®>0. 

lfk = 30® and ?i = ® we obtain Si (30®) = 30®, S2(®) = 30® and it follows 
S2(®)-Si(30®) = 0. 

\ik = 30®2 (we remind that we use here a shorten way to write down this in- 
finite number, the complete record is 30® ^2®°) and « = ® we obtain Si (30®2) = 
30®2, S2(®) = 30® and it follows 

S2 (® ) - Si (30®2) = 30® - 30®2 = -2 < 0. 

Analogously, the expression -^44 can be calculated. 

Let us consider now the famous divergent series with alternate signs S3 = 1 — 
1 + 1 — 1 + ... In literature there exist many approaches giving different answers 
regarding the value of this series (see [5|). All of them use various notions of 
average. However, the notions of sum and average are different. In our approach, 
we do not appeal to average and calculate the required sum directly. To do this we 



should indicate explicitly the number of items, k, in the sum. Then 

0, iik = 2n 



53(.) = 1-1 + 1-1 + 1-1 + 1 ., , ^ if, = 2. + l, 

and it is not important whether k is finite or infinite. 

We conclude this subsection by studying the series ^4 = YX^i ^ converging to 
one. The new approach allows us to give a more precise answer. Due to Postulate 3, 
the formula 

can be used directly for infinite k, too. For example, if ^ = 3® then 

3® 1 1 

where J^ is infinitesimal. Thus, the traditional answer YX=i ^ — ^ ^^ i^^^ ^ finite 
approximation to our more precise result using infinitesimals. The traditional nu- 
meral systems do not allow us to distinguish results of the sums for infinite values 
of k. More examples can be found in |[T3l . 

Thus, if one is able to calculate a partial sum of a series S, he/she can use the 
formula applied for this calculation to evaluate at the Infinity Computer sums S{k) 
with k items for finite and infinite values of k and finite, infinite, and infinitesimal 
values of S{k) and to use the obtained results in further calculations. 

3.2 Computing expressions with infinite and infinitesimal arguments 

In the traditional analysis, the concept of the limit has been introduced in order 
to avoid difficulties that one faces when he/she wants to evaluate an expression 
at infinity or at a point x infinitely close to a point a. If limY^^,/(;c) exists, then 
it gives us a very poor - just one value - information about the behavior of f{x) 
when X tends to a. 

Now we can obtain significantly more rich information using the Infinity Com- 
puter independently on the fact of existence of the limit. We can calculate f{x) 
directly at any finite, infinite, or infinitesimal point expressible in the new posi- 
tional system even if the limit does not exist. Thus, limits can be substituted by 
precise numerals f{a) that are different for different infinite, finite, or infinitesi- 
mal values of X = a. This is very important for practical computations because 
this substitution eliminates indeterminate forms, i.e., again the Infinity Computer 
should not stop its calculations as ti'aditional computers are forced to do when they 
encounter indeterminate forms. 

Example 3.1. In the traditional analysis, the following two limits 

Um (/ + n.5x^ + lOi"") = +00, lim (/ + II.Sjc^) = +00. 



give us the same result, +00, in spite of the fact that for any finite x the difference 
between the two expressions is equal to quite a large number 

The new approach allows us to calculate exact values of both expressions, x^ + 
11. 5x^ + 10^°*^ and x^ + 11.5x^, at any infinite (and infinitesimal) x expressible in 
the chosen numeral system. For instance, the choice of x = 3®^ gives the values 

{3®y + 11.5(3®^)^ + 10'"° = 81®*^103.5®^10i°°®° 

and 81® 103.5® , respectively. Consequently, one obtains 

81®^103.5®'^10^°°®° - 81®*^103.5®'^ = 10^"°. D 

An additional advantage of the usage of the Infinity Computer arises in the fol- 
lowing situations. Suppose that we have a computer procedure calculating f{x), we 
do not know the corresponding analytic formulae for f{x), for a certain argument 
a the value f{a) is not defined (or a traditional computer produces an overflow or 
underflow message), and it is necessary to calculate the limv^a/(x). Traditionally, 
this situation requires a human intervention and an additional theoretical investi- 
gation whereas the Infinity Computer is able to process it automatically working 
numerically with the expressions involved in the procedure. It is sufficient to cal- 
culate f{x), for example, at a point x = a + ®^^ in cases of finite a ov a = and 
X = ® in the case when we are interested in the behavior of /(x) at infinity. 

Example 3.2. Suppose that we have two procedures evaluating /(x) = ^-t^ and 
g{x) = ^. Obviously, /(O) and g(0) are not defined and it is not possible to calcu- 
late limv^o/(-^), liinx^oo/(x) and \mix^Qg{x), \mix^oog{x) using traditional com- 
puters. Then, suppose that we are interested in evaluating the expression 

h{x) = {f{x)-l)-g{x). 

It is easy to see that h{x) = 34 for any finite value of x. On the other hand, the 
following limits 

lim/i(x) = (lim/(x) — 2) • limg(x), 

lim h{x) = (lim /(x) — 2) • lim g{x) 

X-^a° .V— >oo .V— )-oo 

cannot be evaluated. The Infinity Computer can calculate h{x) numerically for 
different infinitesimal and infinite values of x obtaining the same result that takes 
place for finite x. For example, it follows 

/®^ + 2® \ 34 , 

h{®)=[ ^ -2 •g = (® + 2-2)-34®-i=34. D 
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It is worthy to notice that expressions can be calculated by the Infinity Com- 
puter for infinite and infinitesimal arguments, even when their limits do not exist, 
thus giving a very powerful tool for studying divergent processes. 

Example 3.3. The limit lim„^+oo/(«), where f{n) = (— 5)"«, does not exist. How- 
ever, we can calculate the expression (— 5)"?i for different infinite values of n. For 
instance, since it can be easily proved that grossone is even (see |l8][l5l), for n = ® 
it follows /((D) = 5®® and for n = ® - 1 we have /(® - 1) = -5®-\® - 1). □ 

Thus, these new computational possibilities of the Infinity Computer allow one 
both to avoid calculating limits theoretically and to increase the accuracy of nu- 
merical computations. 

3.3 Usage of infinitesimals for solving systems of linear equations 

Very often in computations, an algorithm performing calculations encounters a sit- 
uation where the problem to divide by zero occurs. Then, obviously, this operation 
cannot be executed. If it is known that the problem under consideration has a solu- 
tion, then a number of additional computational steps trying to avoid this division 
is performed. A typical example of this kind is the operation of pivoting used when 
one solves systems of linear equations by an algorithm such as Gauss-Jordan elim- 
ination. Pivoting is the interchanging of rows (or both rows and columns) in order 
to avoid division by zero and to place a particularly 'good' element in the diagonal 
position prior to a particular- operation. 

The following two simple examples give just an idea of a numerical usage of 
infinitesimals and show that the usage of infinitesimals can help to avoid pivoting 
in cases when the pivotal element is equal to zero. We emphasize again that the 
Infinity Computer (see lH) works with infinite and infinitesimal numbers expressed 
in the positional numeral system (|6]l, ^ numerically, not symbolically. 

Example 3.4. Solution to the system 



1' 

2 2 






= 


" 2 " 
2 



is obviously given by ;cj = — 1 , ;c2 = 2. It cannot be found by the method of Gauss 
without pivoting since the first pivotal element a\i =0. 

Since all the elements of the matrix are finite numbers, let us substitute the 
element ayi = by ®^ and perform exact Gauss transformations without pivoting: 
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® 
-2® + 2 
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2® 
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I-® 

1 



1-® 
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It follows immediately that the solution to the initial system is given by the finite 
parts of numbers — 1 + —K^ and 2 



1-® 



1-®- 



We have introduced the number ®^ once and, as a result, we have obtained 
expressions where the maximal power of grossone is one and there are rational 
expressions depending on grossone, as well. It is possible to manage these rational 
expressions in two ways: (i) to execute division in order to obtain its result in the 
form (l6]l, dV]); (ii) without executing division. In the latter case, we just continue 
to work with rational expressions. In the case (i), since we need finite numbers 
as final results, in the result of division it is not necessary to store the parts Cp(if 
with p < —I. These parts can be forgotten because in any way the result of their 
successive multiplication with the numbers of the type ci®' (remind that 1 is the 
maximal exponent present in the matrix under consideration) will give exponents 
less than zero, i.e., numbers with these exponents will be infinitesimals that are not 
interesting for us in this computational context. 

Thus, by using the positional numeral system ^, ^ with the radix grossone 
we obtain 



1 




® 
1 



2® 

-4®+2 
-2®+2 



® 
1 



2® 
2®"+l®" 



2®-® -(2®"!® 
2®°l®"i 



The finite parts of numbers -1®° and 2®°!®" 
provide the required solution. 

Example 3.5. Solution to the system 
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-I®o 
2®°l®"i 



- 1 and 2 respectively, then 
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is the following: ;cj = 2, ;c2 = — 2, and ^3 = 1. The coefficient matrix of this system 
has the first two leading principal minors equal to zero. Consequently, the first 
two pivots, in the Gauss transformations, are zero. We solve the system without 
pivoting by substituting the zero pivot by ®^ ' , when necessary. 
Let us show how the exact computations are executed: 
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-2®2 + 3® 
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4®'-7®+l 



4®-+®+3 J 

It is easy to see that the finite parts of the numbers 

8®2 + 2® 



1 
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8® +2® 
4®'+®+3 
-8®'+10® 
4®>®+3 
4®--7®+l 
4®'+®+3 



4®^ + ® + 3 

-8®^ + 10® _ 

4®2 + ® + 3 ~ 

4®2 - 7® + 1 



4®^ + ® + 3' 



-2 + - 



4®^ + ® + 3 
coincide with the corresponding solution x\ 



1 



12® + 6 

4®V® + 3' 
8® + 2 



4®^ + ® + 3 
'., xl = —2, and;c3 = 1. 
In this procedure we have introduced the number ®^ two times. As a result, 
we have obtained expressions where the maximal power of grossone is equal to 2 
and there are rational expressions depending on grossone, as well. By reasoning 
analogously to Example 13.51 when we execute divisions, in the obtained results it 
is not necessary to store the parts of the type Cp®'^ ,p < —2, because in any way 
the result of their successive multiplication with the numbers of the type C2®^ will 
give finite exponents less than zero. That is, numbers with these exponents will be 
infinitesimals that are not interesting for us in this computational context. Thus, by 
using the positional numeral system ^, (|7]l, we obtain 
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1 





® 
-2®^ - ® 

1 



® 
-2®2 + 3® 

4®'-7®+l 
4®'+®+3 
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® 


® 
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-2®2 - ® 


-2®2 + 3® 





1 


l®"-2®-^ 



^ does not contain the part of the type c_2® ^ 



Note that the number l®''-2® 

because the coefficient c_2 obtained after the executed division is such that c_2 

Then we proceed as follows 
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l®"-2® 



-1 
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1 



l®°-2® 



-1 



The obtained solutions x\ = 2 and X2 = —2 have been obtained exactly without 
infinitesimal parts and JC3 = 1 is derived from the finite part of 1 ®°-2®^ ' . □ 

We conclude this section by emphasizing that zero pivots in the matrix are 
substituted dynamically by ®^ . Thus, the number of the introduced infinitesimals 
®^' depends on the number of zero pivots. 
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4 New computational possibilities for mathematical 
modelling 

The computational capabilities of the Infinity Computer allow one to construct 
new and more powerful mathematical models able to take into account infinite and 
infinitesimal changes of parameters. In this section, the main attention is given 
to infinitesimals that can increase the accuracy of models and computations, in 
general. It is shown that the introduced infinitesimal numerals allow us to formalize 
the concept 'point' and to use it in practical calculations. Examples related to 
computations of probabilities and areas (and volumes) of objects having several 
parts of different dimensions are given. 

4.1 Numerical representations of points at an interval 

We start by reminding traditional definitions of the infinite sequences and sub- 
sequences. An infinite sequence {an],a„ G A,n G N, is a function having as the 
domain the set of natural numbers, N, and as the codomain a set A. A subsequence 
is a sequence from which some of its elements have been removed. 

Theorem 4.1. The number of elements of any infinite sequence is less or equal 

to®. 

Proof. It has been postulated that the set N has (D elements. Thus, due to 
the sequence definition given above, any sequence having N as the domain has ® 
elements. 

The notion of subsequence is introduced as a sequence from which some of its 
elements have been removed. Thus, this definition gives infinite sequences having 
the number of members less than grossone. □ 

It becomes appropriate now to define the complete sequence as an infinite se- 
quence containing ® elements. For example, the sequence of natural numbers is 
complete, the sequences of even and odd natural numbers are not complete. One of 
the immediate consequences of the understanding of this result is that any sequen- 
tial process can have at maximum ® elements and, due to Postulate 1 , it depends 
on the chosen numeral system which numbers among ® members of the process 
we can observe. 

By using the introduced, more precise than the traditional one, definition of 
sequence, we can calculate the number of points of the interval [0, 1), of a line, and 
of the A'^-dimensional space. To do this we need a definition of the temi 'point' 
and mathematical tools to indicate a point. Since this concept is one of the most 
fundamental, it is very difficult to find an adequate definition. If we accept (as is 
usually done in modem Mathematics) that a point A belonging to the interval [0, 1) 
is detemiined by a numeral x, ;c G S, called coordinate of the point A where § is 
a set of numerals, then we can indicate the point A by its coordinate x and we ai^e 
able to execute the required calculations. 
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It is worthwhile to emphasize that we have not postulated that x belongs to the 
set, R, of real numbers as it is usually done, because we can express coordinates 
only by numerals and different choices of numeral systems lead to various sets 
of numerals. This situation is a direct consequence of Postulate 2 and is typical 
for natural sciences where it is well known that instruments influence the result of 
observations. It is similar as to work with a microscope: we decide the level of the 
precision we need and obtain a result which is dependent on the chosen level of 
accuracy. If we need a more precise or a more rough answer, we change the lens 
of our microscope. 

We should decide now which numerals we shall use to express coordinates of 
the points. Different variants can be chosen depending on the precision level we 
want to obtain. For example, if the numbers < ;c < 1 are expressed in the form 
^-,/7 E N, then the smallest positive number we can distinguish is 4^. Therefore, 
the interval [0, 1) contains the following ® points 

1 2 ®-2 (D-l 

' ®' ®' ■■■ ® ' ® ■ 

Then, due to Theorem 14.11 and the definition of sequence, ® intervals of the form 
[a — l,a),a G N, can be distinguished at the ray x > 0. Hence, this ray contains ® 
points and the whole line consists of 2®^ points. 

If we need a higher precision, within each interval 

we can distinguish again ® points and the number of points within each interval 
[a — 1 , a) , a G N, will become equal to ® . Consequently, the number of the points 
of this kind on the line will be equal to 2® . 

Continuing the analogy with the microscope, we can also decide to change 
our microscope with a new one. In our terms this means to change the numeral 
system with another one. For instance, instead of the numerals considered above, 
we choose a positional numeral system to calculate the number of points within the 
interval [0, 1) expressed by numerals 

(.a_ia_2---a-(a)-i)'3-(i))fe- (13) 

Theorem 4.2. The number of elements of the set of numerals f l73l ) is equal to b®. 

Proof. The proof is obvious and is so omitted. D 

Corollary 1. The number of points expressed by numerals 

{a^-iaQ-2 ■ . .aiao.a_ia_2 • ■■<^-{m-i)a-m)b (14) 

is equal to b^®. 

Proof. The corollary is a straightforward consequence of Theorem 14.21 D 
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Figure 1: What is the probability that the rotating disk stops in such a way that the 
point A will be exactly in front of the arrow? 



4.2 Applications in probability theory and calculating volumes 

A formalization of the concept 'point' introduced above allows us to execute more 
accurately computations having relations with this concept. Very often in scien- 
tific computing and engineering it is required to construct mathematical models for 
multi-dimensional objects. Usually this is done by partitioning the modelled object 
in several parts having the same dimension and each of the parts is modelled sep- 
arately. Then additional efforts are made in order to provide a correct functioning 
of a model unifying the obtained sub-models and describing the entire object. 

Another interesting applied ai^ea is linked to stochastic models dealing with 
events having probability equal to zero. In this subsection, we first show that the 
new approach allows us to distinguish the impossible event having the probabil- 
ity equal to zero (i.e., P{0) = 0) and events having an infinitesimal probability. 
Then we show how infinitesimals can be used in calculating volumes of objects 
consisting of parts having different dimensions. 

Let us consider the problem presented in Fig. [T]from the traditional point of 
view of probability theory. We start to rotate a disk having radius r with the point 
A mai^ked at its border and we would like to know the probability P{E) of the 
following event E: the disk stops in such a way that the point A will be exactly in 
front of the arrow fixed at the wall. Since the point A is an entity that has no extent 
it is calculated by considering the following limit 



P{E) 



lim 



0. 



where h is an arc of the circumference containing A and Inr is its length. 

However, the point A can stop in front of the an^ow, i.e., this event is not im- 
possible and its probability should be strictly greater than zero, i.e., P{E) > 0. The 
new approach allows us to calculate this probability numerically. 

First of all, in order to state the experiment more rigorously, it is necessary to 
choose a numeral system to express the points on the circumference. This choice 
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will fix the number of points, K, that we are able to distinguish on the circumfer- 
ence. Definition of the notion point allows us to define elementary events in our 
experiment as follows: the disk has stopped and the anow indicates a point. As a 
consequence, we obtain that the number, N{Q.), of all possible elementary events, 
ei, in our experiment is equal to K where Q. = U,Jj e, is the sample space of our 
experiment. If our disk is well balanced, all elementary events are equiprobable 
and, therefore, have the same probability equal to jwm- Thus, we can calculate 
P{E) directly by subdividing the number, N{E), of favorable elementary events by 
the number, K = N{D.), of all possible events. 

For example, if we use numerals of the type 4r , / G N, then K = ®. The number 
N{E) depends on our decision about how many numerals we want to use to repre- 
sent the point A. If we decide that the point A on the circumference is represented 
by m numerals we obtain 

where the number ^ is infinitesimal if m is finite. Note that this representation is 
very interesting also from the point of view of distinguishing the notions 'point' 
and 'arc'. When m is finite than we deal with a point, when m is infinite we deal 
with an arc. 

In the case we need a higher accuracy, we can choose, for instance, numerals 
of the type /®^^, 1 < / < ®^, for expressing points at the disk. Then it follows 
K = ®^ and, as a result, we obtain P{E) = m®^^ > 0. 

This example with the rotating disk, of course, is a particular instance of the 
general situation taking place in the traditional probability theory where the prob- 
ability that a continuous random variable X attains a given value a is zero, i.e., 
P{X = a) =0. While for a discrete random variable one could say that an event 
with probability zero is impossible, this can not be said in the case of a continuous 
random variable. As we have shown by the example above, in our approach this 
situation does not take place because this probability can be expressed by infinites- 
imals. As a consequence, probabilities of such events can be computed and used in 
numerical models describing the real world (see ||T31 for a detailed discussion on 
the modelling continuity by infinitesimals in the framework of the approach using 
grossone). 

Moreover, the obtained probabilities are not absolute, they depend on the accu- 
racy chosen for the mathematical model describing the experiment. There is again 
a straight analogy with Physics where it is not possible to obtain results that have 
a precision higher than the accuracy of the measurement of the data. We also can- 
not obtained a precision that is higher than the precision of numerals used in the 
mathematical model. 

Let us now consider two examples showing that the new approach allows us to 
calculate areas and volumes of a more general class of objects than the traditional 
one. In Fig.|2]two figures aie shown. The traditional approach tells us that both of 
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Figure 2: It is possible to calculate and to distinguish areas of these two objects 



1 



1 1 

Figure 3: New possibilities for calculating volumes of objects 



them have area equal to one. In the new approach, if we use numerals of the type 
/®^^ / G N, to express points within a unit intei^val, then the unit interval consists of 
® points and in the plane each point has the infinitesimal area ®^' • ®^' = ®^^. As 
a consequence, this value will be our accuracy in calculating areas in this example. 
Suppose now that the vertical line added to the square at the right figure in Fig. [2] 
has the width equal to one point. Then we are able to calculate the area, 52, of the 
right figure and it will be possible to distinguish it from the area. Si , of the square 
on the left 



1-1 = 1, 



S2 = 1 • 1 + 1 • ® 



-1 



1®"1® 



-1 



If the added vertical line has the width equal to three points then it follows 

^2 = l-l + 3-®"^ = 1®*'3®"^ 

The volume of the figure shown in Fig.[3]is calculated analogously: 

V = 1 . M + M . (D-i + 1 . (D-i . ®-i = l®°l®-4®"^ 

If the accuracy of the considered numerals of the type /®^ , ' G N, is not sufficient, 
we can increase it by using, for instance, numerals of the type /®^^, 1 < / < ®^. 
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Then the unit interval consists of ® points and at the plane each point has the 
infinitesimal area ®^ • ®^ = ®^ . As a result, by a complete analogy to the 
previous case we obtain for lines having the width, for instance, equal to five points 
in all three dimensions that 

S2 = I ■ I + 5 -Q-^ = l®°5®-^, 

V = l.\.\ + l.l.5-®-^ + l-5-®-^-5-®-^ = 1(D°5®"225®"1 

Finally, we conclude the paper by a remark that thanks to the new approach 
it becomes possible to measure fractal objects at infinity. A detailed discussion 
on measuring fractals by the introduced infinite and infinitesimal numbers can be 
found in |[T3l [141. 
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